Magnetic induction in a turbulent flow of liquid sodium: mean behaviour and 
slow fluctuations 
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We study the flow response to an externally imposed homogeneous magnetic field in a turbulent 
swirling flow of liquid sodium — the VKS2 experiment in which magnetic Reynolds numbers Ft m 
up to 50 are reached. Induction effects are larger than in the former VKS1 experiment U 0). 
At Rm larger than about 25, the local amplitude of induced field components supersedes that 
of the applied field, and exhibits non-Gaussian fluctuations. Slow dynamical instationarities and 
low-frequency bimodal dynamics are observed in the induction, presumably tracing back to large 
scale fluctuations in the hydrodynamic flow. 
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I. INTRODUCTION 

The induction of a magnetic field by turbulent flows 
of electrically conducting liquids has received a great in- 
terest over the last decade, motivated by a better under- 
standing of astrophysical and geophysical dynamos — i.e. 
the generation of a self-sustained magnetic field by the 
flow of a conducting fluid (|3|). For instance, Marie et at 
(0), Bourgoin et al. (@) and Petrelis et al. Q) study 
the so-called a- and w-effects. Spence et al. (|5|) recently 
evidence induction effects that can only be attributed to 
the temporal fluctuations of the flow. In the same ex- 
periment, Nornberg et al. ((61) evidence periods of inter- 
mittent growth of the induced magnetic field. Volk et al. 
(0) study the long-time fluctuations of induction profiles 
in liquid gallium. Stepanov et al. (Q) study the tur- 
bulent a-effect in a non-stationary flow. Finally, Sisan 
et al. (@) study the influence of the Lorentz force on 
the turbulent dissipation and on the mean flow. A com- 
mon feature of these studies is the use of an apparatus 
in which the instantaneous flow significantly differs from 
its time average. An interesting open issue is to quantify 
the role of fluctuations in the magnetic induction effects. 
Indications can be obtained by comparing induction mea- 
surements for different configurations in the same general 
experimental set-up. 

This is one of the goals of the VKS (von Karman 
sodium) experiment which studies the magnetohydrody- 
namic behavior of the sodium flow generated inside a 
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cylinder by counter rotation of various impellers (fig- 
ure [TJ. In a first series of experiments in 2000-2002 
— VKS1 — , external steady magnetic fields are applied 
either along the cylinder axis or perpendicular to it 
{H; 0; 0). Modifications of this original set-up are mo- 
tivated by the realization that electrical boundary condi- 
tions play an important role in the strength and geometry 
of magnetic induction (fiol ; Hit U3 . fl3l ) and that fluctua- 
tions, as already noted, may have a leading effect. In a 
second evolution of the set-up — VKS2 — , the flow cell 
is enclosed within a layer of sodium at rest and thick 
copper casing resulting in different electrical boundary 
conditions. 



This article presents the flow response to an externally 
imposed homogeneous magnetic field in the VKS2 ex- 
periment. Both VKS1 and VKS2 set-ups are described 
and compared in section [TTJ We report in section IIIII in- 
duction characteristics in the presence of an externally 
applied field and compare them to that of VKS1. We 
identify and study a bimodal dynamics, for which the 
field abruptly changes between two states and stay there 
for variable durations, resulting in a very long timescales 
dynamics. In these regimes, the probability density func- 
tions of some of the components of the induced field, 
particularly the one aligned with the applied field, are 
non-Gaussian. Complementary observations are made 
in section IIVI when a thin annulus is placed in the mid 
shear layer. Finally, we discuss in section IVl some impli- 
cations of our results regarding the dynamo capacity of 
VKS flows. 
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II. THE VKS EXPERIMENTAL SET-UP AND 
CONFIGURATIONS 




FIG. 1 Sketch of the VKS2 experimental set-up. The inner 
and outer cylinders and the optional mid-plane annulus are 
made of copper (in gray). Other parts are stainless-steel. The 
dimension are given in millimeter (left) and normalized by the 
inner cylinder radius R (right). The 3D Hall probe is located 
at point P, 0.25i? from the axis. The angles characterizing 
the orientation of the induced field are displayed on the right. 
Magnetic measurements are made either at point P or at point 
Q when the annulus is present. 

Up to now, all VKS-experiments concern a swirling 
flow produced by two bladed facing impellers distant of 
371 mm in a 412 mm diameter casing. The VKS1 flow 
region is bounded by a 10 mm copper shell and a 20 mm 
steel casing. The VKS2 evolution (figure[T]) basically con- 
sists of putting VKS1 inside a larger sodium volume and 
a 45 mm copper casing, and designing new impellers. 
With respect to VKS1, the motor power has been in- 
creased from 150 kW to 300 kW and the volume of the 
conducting domain is twice larger. The outer radial layer 
of sodium at rest surrounding the flow is of thickness 
OAR. The counter-rotating impellers producing the flow 
are made of stainless steel, are of radius 0.75i? and have 
8 curved blades of 41 mm height and 190 mm curvature 
radius. They are labeled 'TM73' in reference to the study 
of Ravelet et al. (fl3l ) . The velocity unit is based on the 
impellers rotation frequency F and includes the stirring 
efficiency factor V = 0.6 of the TM73 impellers. This 
yields an integral Reynolds number Re = V2ttR 2 F/v. 
The integral kinetic Reynolds number is of the order of 
5 x 10 6 . In addition, the shear-layer instability is a strong 
source of flow instationarity. 

Regarding the magnetic induction, the control param- 
eter of the problem is the magnetic Reynolds number 
R m which compares the stretching of the magnetic field 
by velocity gradients to the magnetic diffusion. We de- 
fine R m = V HQ<j2nR 2 F . The linear relation between the 
rotation frequency and the magnetic Reynolds number 
is R m ~ 1.88F at 120°C, and R m ~ 1.73F at 150°C. 
A temperature regulation has been installed, through oil 
circulating in the outer cylinder, in order to perform long 



time measurements in stationary regimes with temper- 
ature fixed in the range 110 to 160°C (recall that the 
electrical conductivity of sodium varies significantly in 
the neighborhood of its melting temperature). The max- 
imum magnetic Reynolds number is of the order of 50, 
compared to about 35 in VKS1. 

In order to study the MHD response of the flow, we 
apply a transverse field Bo = Boe y with a pair of coils. 
They are not in Helmholtz configuration and the inho- 
mogeneity of Bo, defined on the flow volume, can reach 
up to 20% at the outer flow boundary. We measure the 
three components of the magnetic field B with a 3D-Hall 
probe set in the equatorial plane inside the flow, 50 mm 
from the axis, i.e. at point P in figure [TJ In this case, 
the closest VKS1 configuration |l|) has impellers TM70 
of radius 0.75R and baffles on the cylinder wall. The ap- 
plied field is less than 3 gauss, too weak a field to modify 
the flow: we have checked that the induced field varies 
linearly with the applied field. 



III. RESPONSE TO AN EXTERNALLY APPLIED FIELD 

Under the imposed field Bq = 2.7 G, we study the 
dynamics of the three components of the dimensionless 
induced field b = (B— B )/-Bo and of its orientation at P. 
We define 9 a latitude with respect to the x — y plane, 
and 4> the longitude in this plane, with respect to the 
applied field (y axis) so that b = ||b|| [cos# (— sincAe x + 
cos0e y ) + sin#e z ]. With these conventions, the applied 
field at P is of unit norm with orientation 9 = 0, cj> = 0. 

We first present our observations about the mean mag- 
netic induction properties of the flow and thereafter dis- 
cuss the temporal dynamics of the signal. 



A. Mean induction 

We plot in figure [2] the evolution of the mean (a) and 
standard deviation (b) of the components of b. The mean 
vertical field (b z ) is zero, which can be explained by the 
symmetries of the time-averaged flow and of the applied 
field SEi). The induced field is dominated by the axial 
component (b x ) and once R m > 20, the transverse com- 
ponent (by) becomes greater than the applied field. Due 
to limitations of the mechanical seals, we are limited to 
F > 6Hz (R m > 11). 

Results for the mean induction can be compared to the 
VKS1 experiment {J @). Recall that the main changes 
concern the electrical boundary condition due the added 
radial blanket of sodium, and the modification of the 
impellers. For R m > 12, we observe that the behav- 
ior of (h)(R m ) is linear, as in VKS1, but with a larger 
slope. We find (b x )/B = R m /R* m with l/R* m = 1/13, 
a value 25% greater than in VKS1. We also observe 
that the induced field no longer saturates for R m > 25: 
the induction is still growing at high R m , leading to 
(b x ) ~ 3.3 ; (by) ~ 2.3 at the highest R m ~ 50, 
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FIG. 2 Mean value (a) and rms amplitude (b) of the induced 
magnetic field vs. R m , for the reduced magnetic field b = 
(B - Bo) /So induced at point P. (a): b x , (□): by, (A): b z . 
and (+): ||b||. 



while the corresponding VKS1 experiment saturates at 
(6a,) ~ 1.5 ; (by) ~ 0.4 for R m > 25. Actually, for 
Rm ^ 25, the field component induced in the direction 
of the applied field is on average larger than the applied 
field although the system shows no sign of dynamo-self 
generation — the induced field vanished as the applied 
field is turned off, save for induction due to the local 
Earth magnetic field. 

Finally, we observe that the standard deviations of the 
three components also increase linearly with R m . Again 
there is no indication of a saturation for the larger R m 
values. The rms amplitude are of the order of 50% of 
the mean induced field. We note that for R m > 45, fluc- 
tuations in each component of the induced field exceed 
the magnitude of the applied field Bo- In addition, the 
fluctuations of the axial component are twice larger than 
the fluctuations of the other components (figure In 
comparable VKS1 measurements the amplitude of fluc- 
tuations are the same for all components, with value 
a(h) ~ 0.5 at R m ~ 25. 



B. Temporal and spectral dynamics 

In order to describe the time dynamics of the magnetic 
induction, we focus on measurements at R m = 40 (F — 
22 Hz). We plot in figure [3] the time series, probability 
density function (PDF) and spectra of the components of 
the induced field as well as the 9 and <j> field orientation 
angles. As could be inferred from the magnitude of the 
standard deviations, one observes strong fluctuations in 
the induced field. 



Whereas in VKS1 all components of the magnetic 
field exhibit almost Gaussian PDF, we now observe non- 
Gaussian statistics for two components: b x and b y (fig- 
ure [3p-g). The component aligned with the applied field 
by (figure [3Jd and [3g) exhibits an exponential tail of the 
form H(b y ) oc exp(— 1.39 x by). Its centered and reduced 
PDF does not depend on R m . Upon closer inspection of 
figure [3^, one also detects bimodality in the time evo- 
lution of the b x component. Relatively long periods of 
high induced field (e.g. around t = 600 and 800F -1 ) are 
followed by periods when b x is around zero {e.g. around 
t = 700 and 92QF~ 1 ). This bimodality is analyzed in 
section lHI.CI 

We also show in figure [3l the time spectrum of the 
three components of b. At high frequency, the fall-off 
of the power spectra is steeper than the —11/3 power 
law decrease expected for the magnetic dissipative range 
in a Kolmogorov- like turbulence (|15l ). i.e. for frequencies 
/ Z Rm 3/4 -F ~ 20.F at R m = 40. This steeper slope was 
also observed in previous VKS1 measurements, although 
the expected Kolmogorov behavior was observed in von 
Karman gallium flows stirred by rugose disks (jig). 

At low frequencies, roughly between F/60 < / < F, 
the spectrum behaves approximately as a power law with 
an exponent —1 for b x . For the two other components, 
the exponent is of the order of —0.5. This type of spec- 
tral behavior is indicative of long-time evolutions in the 
magnetic induction. We believe that it traces back to 
slow changes in the underlying hydrodynamic flow, for in- 
stance the chaotic dynamics of the azimuthal shear layer 
observed in water prototype flows (fi"7l : HH: liH). It is also 
observed for induction measurements in von Karman gal- 
lium flows (|7|) and is already present in VKS1 measure- 
ments. 



C. Bimodal analysis: High-State and Low-State 

We consider the bimodal dynamics illustrated above on 
the axial component of the magnetic field. These jumps 
between two states, which arise irregularly in time, ap- 
pear even clearer when one plots the time evolution of the 
orientation of the induced field b, as in figure [3Ji-e. In 
the periods during which the b x component is large, the 
field orientation has smaller fluctuations whereas during 
periods for which b x is almost null, the orientation has 
strong fluctuations. Hence, the PDF of the axial induced 
component exhibit a clear bimodal behavior shown in fig- 
ure |3f and the PDFs of the field orientation are peaked 
but with large wings (figure Oi) : the PDF of is sym- 
metric around the mean value 0, while the PDF of 4> has 
a peak around —60°, close to the mean orientation, but 
is skewed towards positive angles. 

In the following, we separate the time signals into two 
parts with a threshold criterion 6_ on b x , illustrated in 
figure Hi for F = 22Hz (R m = 40) with = 1.9. We 
define the L— state (resp. H— state) as the subset when 
b x < 6_ (resp. b x > 6_) during a time larger than t m i n ~ 
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FIG. 3 (a-e) Temporal signal of the components and of the orientations of b at R m = 40. (f,g) Probability density function 
(PDF) of b x and b y at R m = 40 (solid line). The dash-dotted line corresponds to the PDF in the H— states and the dashed 
line to the PDF in the L— states — see section llll.Cl for details about the separation. Solid black lines are Gaussians. (h) Same 
plots for orientation angles 9 and <j>. (i) Time spectrum of b x (blue), b y (red) and b z (green) at R m = 40. 



F _1 . The remaining pieces are taken as L— state (resp. 
H— state) if the previous and the following signals are in 
the L— state (resp. H— state). This allows short magnetic 
field fluctuations below 6_ to belong to the H— state and 
vice versa. We also build a binary state signal with value 
1 in the H— states, and in the L— states (figure 0^,). 

We observe in figure [3f that the PDF of b x in each 
state (dashed-dotted and dashed lines) are close to Gaus- 
sian fits of mean 0.66 and standard deviation 1.21 in the 
L— state with corresponding values equal to 3.21 and 1.59 
in the H— state (at R m — 40). Note that the shape of 
the PDF of by (figure [3g) and b z (not shown) are not 
affected by the bimodal decomposition. The transverse 
component b y has its most probable value slightly shifted 
to a lower value in the L— state, with again exponential 
tails in both states — a robust feature not affected by the 
bimodal decomposition. The PDF of the angles for the 
two separated states is plotted in figure [3h- The field 
orientations are much more fluctuating in the L— state 
(standard deviation cr(<fr) = 34°) than in the H— state 
(cr(</>) = 14°). The wide tail towards high values of <f> 
belongs to the L— state. 



The bimodal analysis gives some insight into the slow 
evolution of the induction processes in the flow. It is best 
seen when one considers the autocorrelation functions of 
the induced components, i.e., another way to analyze the 
low frequency part of the power spectra. We plot the au- 
tocorrelation functions for b x , b y , b z in figure 0Jd. They 
decrease exponentially at large time lags, and we perform 
a fit of these curves with a sum of exponentials. For 
by (b z gives approximately the same results), the func- 
tion is well fitted by two exponentials of decay times 2.4 
and 0.18F -1 . The b x component behaves differently and 
three exponentials are necessary to fit the curve. The 
corresponding time scales are 9.0, 1.1 and 0.1 If _1 . The 
two last times are of the same order of magnitude as for 
by and b z : one is of the order of the disks period F^ 1 and 
the shortest one is of the order of the period of a blade 
(since there are 8 blades, the period is 0.125F" 1 ). How- 
ever, b x shows a slow evolution with a time scale of the 
order of 9F _1 . This time arises from the bimodal tran- 
sitions. Indeed, when we compute the autocorrelation 
for the binary state signal, we find that two exponentials 
are sufficient for the fit. Their characteristic times are 
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FIG. 4 (a) Bimodal decomposition of b x component: 
H— state (red) and L— state (blue). Binary state signal 
(black), (b) Autocorrelation functions in lin-log scale for the 
induced field components b x (blue), b y (red) and b z (green), 
with nonlinear fits by a sum of exponentials. Black solid line: 
autocorrelation function of the binary state signal. Measure- 
ments for Rm — 40; see text for details. 



9.6 and 1.3 F~ l . We thus recover the two first times of 
b x . The long time-scale, visible only for the binary state 
signal and for the complete b x signals, corresponds to 10 
times the impellers period and gives a bound correspond- 
ing to / = 1/(2tt x 10F" 1 ) ~ F/60 in the spectrum. We 
have checked for various F that it seems to be a hydrody- 
namical time, of the order of magnitude of the coherence 
time scale of the largest coherent structures in the flow, 
created by the shear layer instability (fl7l; Hit). 

Finally, we show in figure [5] the evolutions with the 
magnetic Reynolds number R m of the mean and standard 
deviations for the induced magnetic field in the High- and 
Low-states. In each state, the mean values evolve linearly 
with R m but the components behave in a different way: 
in the H— state (figure , they are close to the mean of 
the total field, with (b x ) > (b y ), whereas in the L— state 
(by) > (b x ) (figured). The fluctuations of the induced 
field (figure [5fc-d) are of the same order of magnitude in 
both states and are isotropic in the L— state. 



IV. INDUCTION WITH AN ADDITIONAL ANNULUS IN 
THE SHEAR-LAYER. 

Extensive testing and visualizations in water prototype 
experiments (fl8r ) have shown that a way to stabilize the 
shear layer in the center of the flow is to introduce an an- 
nulus in the mid plane. Its leading effect is to reduce the 
low-frequency instationarities in the velocity field while 
the smaller scale fluctuations remain unchanged. As an 



FIG. 5 Evolution with Rm for the H and L-states. (a,b) 
mean value in the L-state and in the H state; (c,d) same for 
the standard deviation, (o): b x , (□): b y , (A): b z . and (+): 
llbll. 



illustration, the main structure of the free shear layer con- 
sists of three big fluctuating vortices. In the presence of 
an annulus, each vortex splits into a pair of smaller vor- 
tices that remain mostly attached to the leading edges 
of the annulus, one on each face. We made a series of 
sodium experiments for an annulus with inner diameter 
175 mm inserted along the inner cylinder in the mid- 
plane between the disks. Magnetic measurements are 
then made at point Q, flush with the outer flow wall, 
just behind the annulus (figure[l}. The distance between 
points P and Q is of the order of several magnetic dif- 
fusion length, preventing complete quantitative compar- 
isons. However, several interesting features emerge. 

We show time series and time spectra of the induced 
magnetic field components for two R rn in figure [5] First, 
the time spectra of the b x and b y components level off 
for / < 0.1F indicating that several long-time dynam- 
ical features may have been suppressed, especially at 
low R m . Instead, time-series and spectra do reveal a 
low-frequency oscillation between F/10 and F/5. This 
oscillation dominates the signal at low R m (figure [5^,- 
d). At higher R m (figure [(J-h), this oscillation ap- 
pears combined with lower frequency drifts — especially 
on b z . However, its dynamics is never comparable to the 
characteristic very-low-frequency binary jumps between 
L— state and H— state without annulus : the fluctuations 
are indeed Gaussian — in the sense that the PDF of all 
components bi(t) follow a Gaussian distribution. 

We observe that the evolution of the mean induction 
also differs (figure [7]). The mean amplitude of b still 
depends linearly on R m . It is dominated by the b y com- 
ponent — parallel to Bo- However, the evolution of b z 
and of the standard deviation of b y and b z shows a re- 
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FIG. 6 Temporal signals of the three components of the induced magnetic field and power spectra, with the annulus in the 
midplane. Measurements are done at point Q — see figure Q] b x (blue), b y (red) and b z (green), (a-d): R m = 15. (e-h): 
R m = 42. The dimensional fields are divided by the applied field measured at Q, i.e., Bq = 1.8 G. 
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FIG. 7 Mean value (a) and rms (b) of the induced magnetic 
field vs. R m , with the annulus in the midplane. (o): b x , (□): 
b v , (A): b z and (+): j|b||. 



markable new feature compared to the previous induction 
measurements (including those performed in the VKS1 
configuration): they are no longer linear function of R m . 
Above R m ~ 30 they start growing faster, together with 
the apparition of the very-low-frequency part of the b z - 
spectrum, i.e. below F/10. 



V. DISCUSSION AND CONCLUDING REMARKS 

Comparison with VKS1 experiment. The measure- 
ments performed show that the addition of an surround- 
ing layer of sodium and the modification of the driving 
impellers have lead to an increased induction efficiency as 
compared to former VKS1 runs. Regarding the bound- 
ary condition, this is in agreement with general consider- 
ations and numerical simulations: the effective magnetic 
Reynolds number is increased if currents can develop over 
a larger volume than the velocity domain (fill: fl4l [20: [2lh . 
Regarding the flow generation, it shows the sensitivity 
of the induction to the precise geometry of the velocity 
gra dients linked to the design of the driving impellers 

EiLil- 

Effects of small appendices on temporal fluctuations 
in confined turbulent flows. Another finding is that the 
long-time dynamics observed here, as in many other tur- 
bulent flows in confined geometries, can be significantly 
altered by relatively small mechanical appendices, such 
as the thin annulus we placed in the equatorial plane 
(section H*V|) . We indeed observe that, whereas the mean 
induction remains of the same order of magnitude, the 
slowest time-scale of the induced magnetic field is re- 
duced by one order of magnitude with respect to the 
free case. This is in agreement with the hypothesis that 
such slow evolutions of the magnetic field are linked to 
large scale non-stationarities in the velocity field. In 
von Karman flows, they originate mainly in the dynam- 
ics of the shear layer created by the counter-rotation of 
the driving disks, as evidenced in water prototype mea- 
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surements (flU Il9l ) and induction measurements in liquid 
gallium — for which the very low value of the magnetic 
Reynolds number result in the magnetic induction actu- 
ally behaving as an image of velocity gradients (0). 

Concluding remarks. Some of the observations de- 
scribed above may have implications regarding the self- 
generation of a magnetic field in a von Karman flows. Dy- 
namo predictions are often obtained from kinematic sim- 
ulations in which the velocity field is time-independent, 
fixed to its time average value (v) (r) . In our case (v) (r ) is 
measured from water model experiments. The threshold 
for which dynamo action is predicted varies significantly 
with prescribed boundary conditions, and, to a lesser ex- 
tent, with the numerical procedure. Values range from 
R c m = 46 for (finite axial boundary condition / 5 mm 
copper shell separating the flow and the static conducting 
layer /copper container / no fluid behind the impellers) 
to R c m > 100 for the same previous conditions but with 
fluid in motion behind the impellers (fl3t I22I : [23h . Nu- 
merical studies - kinematic (fill ) and dynamic (|24[ ) - re- 
veal that non linear trends should be detected when the 
magnetic Reynolds number reaches about 70% of the dy- 
namo threshold value. If this can directly be applied to 
the experiment, our results indicate that in the absence 
of the annulus, the critical magnetic Reynolds number, 
if it exists, should be larger than about 70. Note also 
that although we did observe that the induced magnetic 
field at the measurement location can greatly exceed the 
amplitude of the applied field (cf. figure 2), dynamo self- 
generation did not result. 

One related observation is that in order for dynamo 
generation to occur, several induction processes must co- 
operate. For instance, it has been shown that large scale 
phase fluctuations of the eddies of the G. O. Roberts flow 
(|25r). or adding random large scale noise to the Taylor- 
Green flow (|26r ) could significantly increase the dynamo 
threshold (although natural Navier-Stokes fluctuations 
may actually decrease it compared to the kinematic value 
computed from the mean flow |27|). The introduction of 
an inner annulus, known to stabilize large scale fluctua- 
tions in the velocity field, did not lead to dynamo genera- 
tion, but the measurements have shown a clear non- linear 
increase of the fluctuations of induction. 

It remains for future study to determine how crucial 
each of these factors is. Dynamo generation has been 
achieved in the VKS2 experiment (|28T ) for a rather low 
magnetic Reynolds number, around R m =31. The me- 
chanical dynamo configuration corresponds to the second 
configuration (section HV|) of the present article, i.e. with 
an annulus in the mid-plane, but the steel impellers used 
here are replaced by pure iron ones. Since these are of 
identical shape, the hydrodynamic flow is presumably the 
same (we checked that power consumption is identical 
in both flow below dynamo threshold) but the boundary 
conditions are quite different: the magnetic field does not 
penetrate the region behind the disks and magnetic field 
lines have to attach nearly perpendicular to the surface 
of the iron impellers. 
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